namePool = [0.87];
f = [0.87] * 1e12;



[Ip, Im, IT, IQ, IM] = FC_RadiativePower_U('.\SC_resonance_point\SC_',namePool, f);
YW_Plot_Fast_U({f / (1e12), Ip.Ip, f / (1e12), Im.Im, f / (1e12), IT.IT, f / (1e12), IQ.IQ, f / (1e12), IM.IM}, {'-o', '-d', '-v', '-s', '-^b'}, ["\fontsize{28}\fontname{Helvetica}\itP", '\fontsize{28}\fontname{Helvetica}\itM', '\fontsize{28}\fontname{Helvetica}\itT', '\fontsize{28}\fontname{Helvetica}\itQ_e', '\fontsize{28}\fontname{Helvetica}\itQ_m'], ["Frequency (THz)", 'Radiation   power(a.u.)']);

YW_Plot_Fast_U({f / (1e12), IT.ITx, f / (1e12), IT.ITy, f / (1e12), IT.ITz, f / (1e12), IT.IT}, {'-o', '-d', '-v', '-s'}, ["\fontname{Helvetica}\itTx", '\fontname{Helvetica}\itTy', '\fontname{Helvetica}\itTz', '\fontname{Helvetica}\itT'], ...
    ["Frequency (THz)", 'Radiation   power (a.u.)']);
set(gca, 'yscale', 'log');
